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C\l ; ABSTRACT 

We present one-zone chemical evolution models for two dwarf starburst galaxies, NGC 1705 
and NGC 1569. Though especially designed for the inner ~ 1 kpc region, where numerous 
H II regions and most of the stars are observed, the models also account for the presence of 
extended gaseous and dark matter haloes, and properly compute the binding energy of the gas 
heated by supernova explosions. Using information about the past star formation history and 
initial mass function of the systems previously obtained from Hubble Space Telescope opti- 
cal and near-infrared colour-magnitude diagrams, we identify possible scenarios of chemical 
(*f) ' enrichment and development of galactic winds. We assume that the galactic winds are propor- 

tional to the Type II and Type la supernova rates. As a consequence, they do not necessarily 
go to zero when the star formation stops. In order not to overestimate the current metallicity of 
the interstellar gas inferred from H II region spectroscopy, we suggest that the winds efficiently 
remove from the galaxies the metal-rich ejecta of dying stars. Conversely, requiring the final 
mass of neutral gas to match the value inferred from 21 -cm observations implies a relatively 
low efficiency of interstellar medium entrainment in the outflow, thus confirming previous 
findings that the winds driving the evolution of typical starbursts are differential. These con- 
Qh elusions could be different only if the galaxies accrete huge fractions of unprocessed gas at 

late times. By assuming standard stellar yields we obtain a good fit to the observed nitrogen 
^ ■ to oxygen ratio of NGC 1569, while the mean N/O ratio in NGC 1705 is overestimated by the 

^ | models. Reducing the extent of hot bottom burning in low-metallicity intermediate-mass stars 

does not suffice to solve the problem. Localized self-pollution from stars more massive than 
60 M Q in NGC 1 705 and/or funneling of larger fractions of nitrogen through its winds are then 
left to explain the discrepancy between model predictions and observations. Inspection of the 
. log(N/0) vs. log(0/H)+12 diagram for a large sample of dwarf irregular and blue compact 

' dwarf galaxies in the literature favours the latter hypothesis. Yet, the physical mechanisms 

responsible for such a selective loss of metals remain unclear. 

Key words: galaxies: abundances - galaxies: evolution - galaxies: formation - galaxies: 
individual: NGC 1705, NGC 1569 - galaxies: irregular - galaxies: starburst. 
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1 INTRODUCTION 

Low-luminosity galaxies are the most common type of galaxies in 
the nearby universe (Hodge 1971). Dwarf galaxies present a well 
pronounced dichotomy between the two main classes of early-type 
and late-type dwarfs (dwarf irregulars, DIGs, and blue compact 
dwarfs, BCDs). Indeed, most dwarfs are either gas-poor spheroids 
with dominant intermediate-age or old stellar populations, or gas- 
rich, star-forming systems (van den Bergh 1977). Only few tran- 
sition objects are known (Sandage & Hoffman 1991). Here we 
concentrate on DIGs and BCDs. BCDs are currently undergoing 
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vigorous star formation (SF) activity. They are objects of compact 
appearance, with centrally concentrated starburst, gas and star dis- 
tributions. Part of the extended neutral gas may be kinematically 
decoupled from the galaxies (van Zee, Skillman & Salzer 1998). 
DIGs are instead dominated by scattered bright H II regions in the 
optical, while in HI they show a complicated fractal-like pattern of 
shells, filaments and clumps. Typical H I masses are < 10 9 M Q . 

Studying gas-rich starbursting objects helps to better under- 
stand the process of star formation from low-metallicity gas clouds 
and the derivation of the primordial helium abundance with the 
need for a minimum extrapolation to early conditions (e.g., Peim- 
bert & Torres-Peimbert 1974; Izotov & Thuan 1998). Moreover, 
by studying dwarf galaxies, one can complement the information 
available for solar composition stars and get useful insights into the 
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metal dependence of the stellar yields and the winds in Wolf-Rayet 
(W-R) stars. 

Late-type dwarfs have low metallicities, large gas content and 
mostly young stellar populations: all these features indicate they 
are poorly evolved objects and are consistent with two possible ex- 
planations for their nature: (i) either they are newly formed galaxies 
or (ii) they have evolved slowly over the Hubble time. In their pi- 
oneering work, Searle, Sargent & Bagnuolo (1973) discarded the 
youth hypothesis and concluded that extremely blue galaxies must 
have undergone brief but intense bursts of SF separated by long 
quiescent periods (bursting SF). Recent detections of old underly- 
ing stellar populations in most BCDs seem to corroborate Searle's 
et al. (1973) earlier suggestion and reveal at least another burst of 
SF besides the present one, even in the case of the most metal-poor 
BCD known, IZw 18 (Aloisi, Tosi & Greggio 1999; Ostlin 2000; 
Izotov & Thuan 2004). 

Besides the bursting SF mode, gasping (Tosi et al. 1991) or 
mild continuous (Carigi, Colin & Peimbert 1999; Legrand et al. 
2000) SF regimes have been proposed for DIGs and BCDs. The 
gasping regime, in which the interburst periods last significantly 
less than the active phases, probably is the more realistic picture 
for many of them (Schulte-Ladbeck et al. 2001). 

Larson & Tinsley (1978) postulated the existence of a close 
connection between violent dynamical phenomena and bursts of 
SF. Observational evidence of intense ongoing SF in isolated BCDs 
(Campos- Aguilar & Moles 1991) contradicts this hypothesis. How- 
ever, at least some 'isolated' BCDs might have low surface bright- 
ness or pure HI companions (Mendez & Esteban 2000; Noeske 
et al. 2001; Pustilnik et al. 2001; Hoffman et al. 2003). In addi- 
tion to tidal interactions, low-mass stellar or gaseous companions 
might trigger and fuel starbursts by infalling on to the DIG (Koppen 
& Edmunds 1999; cfr. also the 'minor merger' phenomenon de- 
scribed by Sancisi 1999). Another explanation for the intermittent 
mode of SF in DIGs and BCDs was given by Gerola, Seiden & 
Schulman (1980), who applyed the stochastic self-propagating SF 
(SSPSF) model to systems of small size. Yet, the high dispersion of 
the properties of these galaxies cannot be explained in terms of the 
SSPSF mechanism alone (Matteucci & Chiosi 1983). In particu- 
lar, galactic winds (GWs) of variable intensity (Matteucci & Chiosi 
1983; Matteucci & Tosi 1985) releasing metals in different propor- 
tions (enriched winds - Pilyugin 1993; Marconi, Matteucci & Tosi 
1994; Carigi et al. 1995) were introduced in order to reproduce the 
high AY/AO ratio and the distribution of N/O vs. O/H observed in 
extragalactic H II regions. 

Observational evidence for starburst-driven large-scale out- 
flows was established for a number of starburst or young post- 
starburst systems, including IZw 18 (Martin 1996), NGC 1569 
(Waller 1991; Heckman et al. 1995), and NGC 1705 (Meurer et 
al. 1992; Marlowe et al. 1995). This evidence is now growing fast, 
since diffuse and filamentary gas associated with large-scale out- 
flows are detected for an increasing number of starbursting dwarfs 
(Veilleux et al. 2003; Cannon et al. 2004). However, whether the 
final fate of this gas is to eventually escape from the galaxy or to 
cool and recollapse forming the next generations of stars is not yet 
clear. Far Ultraviolet Spectroscopic Explorer (FUSE) observations 
of the O VI A 1032 absorption line arising in the galactic outflow of 
NGC 1705 suggest that the superbubble has begun to blow out of 
the interstellar medium (ISM) in the disk of NGC 1705 (Heckman 
et al. 2001). This superbubble is likely to become a superwind and 
vent its metals and kinetic energy. The bulk of the global ISM (H I) 
in NGC 1705 will be retained (Mac Low & Ferrara 1999). A simi- 
lar scenario should be suitable also for IZw 18 (Martin 1996). Re- 



cently, we have also been provided with the first direct evidence for 
metal-enriched winds from dwarf starburst galaxies: deep Chandra 
imaging of NGC 1569 has led to an estimated mass of oxygen in the 
hot wind similar to the oxygen yield of the current starburst (Mar- 
tin, Kobulnicky & Heckman 2002). The most likely interpretation 
is that the wind carries nearly all the metals expelled by supernovae 
(SNe). A similar scenario emerges also from chemodynamical sim- 
ulations (e.g. Recchi et al. 2004, 2005). 

In this paper we study the chemical evolution of DIGs and 
BCDs. In particular, we compare predictions from models with 
structural parameters applicable to NGC 1705 and NGC 1569 with 
their observed properties. Our choice is motivated by the fact that 
the recent star formation history (SFH), initial mass function (IMF) 
and chemical enrichment of these two galaxies are rather well con- 
strained by the observations. Therefore, we can construct mod- 
els without making arbitrary assumptions about the SFH and IMF 
(which are usually free parameters of the chemical evolution mod- 
els), and assume the values inferred from the observations. As a 
consequence, we are able to better constrain other poorly known 
physical processes which drive galaxy evolution, such as, for in- 
stance, the stellar feedback. Both NGC 1705 and NGC 1569 have 
strong galactic winds, likely triggered by the recent starbursts, and 
we aim at better understanding how these outflows develop. 

The layout of the paper is the following. Observations are pre- 
sented in § 2, the models are described in § 3, model results are 
given in § 4 and discussed in § 5. Finally, the conclusions are drawn 
in § 6. 



2 OBSERVATIONS 

2.1 General properties of dwarf irregular and blue compact 
dwarf galaxies 

2.1.1 Chemical abundances 

Metallicities in DIGs and BCDs are usually derived from the ion- 
ized gas in H II regions. From the optical, O, N, S, Ne, Ar, He, and 
even fainter lines such as those of Fe are measured, while C and Si 
abundances are obtained from the ultraviolet (UV). Oxygen is the 
element with the most reliable determinations, since all its most 
important ionization stages can be observed. The intrinsic uncer- 
tainty is of the order of ~ 0. 1 dex, which rises at 0.2 dex or more 
when the electron temperature cannot be directly determined (Pagel 
1997). HII regions are ionized by newly born massive stars, hence 
showing the metallicity of the ISM at the present time. A possi- 
ble problem with metal-poor H II regions could be self-pollution of 
fresh metals by winds from young massive stars, so that the abun- 
dances inferred from the nebular emission lines would not be rep- 
resentative of that in the local ISM (Kunth & Sargent 1986; Pagel, 
Terlevich & Melnick 1986). However, if the elements ejected from 
massive stars are in high ionization stages, they are not expected to 
significantly contribute to the element abundances as derived from 
optical lines (Pantelaki & Clayton 1987). Observational evidence 
for this comes from the chemical homogeneity of most star-forming 
galaxies, in which the expected spatial variations in oxygen abun- 
dances are not observed, despite the presence of multiple massive 
star clusters (Kobulnicky & Skillman 1997, 1998). 

Recently, a step forward in understanding the physical and 
chemical evolution of dwarf star-forming galaxies has been accom- 
plished thanks to FUSE, which makes possible the determination of 
the metal abundances of the neutral ISM of these galaxies. Prelim- 
inary data analysis points to an offset in metal content between the 
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neutral ISM and H II regions, at least for the galaxies surveyed up 
to now (Aloisi et al. 2005). We naturally find a qualitative explana- 
tion for this offset in the framework of our model (see discussion 
in Sect. 4). However, many uncertainties still affect both the inter- 
pretation of the data and the physics of the model, thus making the 
accurate quantification of the offset a hardly achievable goal. 

Only for the closest galaxies, there is the possibility of trac- 
ing the history of the chemical enrichment of the ISM through 
spectroscopy of the resolved stellar populations (Hill 2004; Tol- 
stoy et al. 2004, and references therein). Photometric estimates of 
the iron stellar abundances can be obtained from colour-magnitude 
diagrams (CMDs), but with much lower reliability. 



2.1.2 Star formation history and initial mass function 

The synthetic CMD method (Tosi et al. 1991; see also Gallart et 
al. 1994; Tolstoy & Saha 1996; Aparicio & Gallart 2004 and refer- 
ences therein) allows one to estimate the epochs and intensities of 
the SF activity in a galaxy, as well as the IMF slope. 

In galaxies outside the Local Group, crowding and magnitude 
limits make it increasingly difficult to resolve the fainter stars and 
shorten accordingly the reachable look-back time. Nevertheless, 
their SFHs have been inferred up to a few Gyr ago (Tosi 2003, 
and references therein). 

The following sections, Sects. 2.2.1 and 2.2.2, are mostly de- 
voted to a detailed description of the SFHs, IMF slopes, and chem- 
ical abundances inferred for NGC 1705 and NGC 1569. In spite of 
being differently classified (as a BCD NGC 1705 and as a DIG 
NGC 1569), they are otherwise similar objects, showing both evi- 
dence of a complex structure and the presence of galactic outflows. 



2.2 Individual objects 

2.2.1 NGC 1705 

NGC 1705 is a nearby, isolated BCD. Its nuclear region contains an 
extremely luminous super star cluster (SSC), responsible for about 
half of the total UV light from the galaxy (Meurer et al. 1992). 
SF is found to occur now, mainly concentrated in the high surface 
brightness part of the galaxy (Annibali et al. 2003). In this region, 
the presence of W-R emission lines suggests ongoing SF during 
the past 5 Myr (Meurer et al. 1992). In Ha, NGC 1705 shows a 
bipolar morphology, suggesting an outflow powered by the recent 
SF activity in the nucleus (Meurer et al. 1992). UV absorption- 
line kinematics bears witness to this too (Heckman & Leitherer 
1997). The gas metallicity, as derived from UV, optical and near in- 
frared (NIR) spectra of H II regions is 12+log(0/H) = 8.36 (Storchi- 
Bergmann, Calzetti & Kinney 1994), corresponding to Z ~ 0.004 
(for Z Q = 0.017), close to the metallicity of the SMC. A higher 
value, 12+log(0/H) ~ 8.46, is reported by Meurer et al. (1992), 
while Heckman et al. (1998) give 12+log(0/H) ~ 8.0. Recently, 
Lee & Skillman (2004) have derived a mean oxygen abundance 
of 12+log(0/H) = 8.21 ± 0.05. They have detected [OIII] A 4363 
in five HII regions of NGC 1705, which makes their oxygen de- 
termination highly reliable, thanks to the direct measurements of 
electron temperatures. The estimated log(N/0) = —1.75 ± 0.06 is 
among the lowest values ever observed for dwarf irregulars, but it 
raises up to —1.63 ± 0.07 if the anomalously low nitrogen abun- 



Table 1. NGC 1705 identity card. 



Quantity 


Observed value 


References 


D 


5.1 ± 0.6 Mpc 


1 


Rh 


1.7 kpc a 


2 


M gSlB 


1.7 x 10 8 


M a 


3 


A^stars 


1.7 x 10 8 


M a 


3 


z 


0.004 


4 




H ii regions 


Neutral ISM 




log(0/H) + 12 


8.46 




3 




8.36 




4 




8.0 




5 




8.21 ±0.05 


7.50±0.01 b 


6,7 


log(N/H) + 12 


6.46 ± 0.08 


5.68 ±0.06 


6,7 


log(Fe/H) + 12 




6.21 ±0.03 


7 


log(N/0) 


-1.75±0.06 c 


-1.82 ±0.07 


6,7 




-1.62 




4 


log(Ne/Q) 


-0.426 ±0.033 




6 



1- Tosi et al. 2001 ; 2- Meurer et al. 1998; 3- Meurer et al. 1992; 4- Storchi- 
Bergmann et al. 1994; 5- Heckman et al. 1998; 6- Lee & Skillman 2004; 
7- Aloisi et al. 2005. 

a Radius and masses were modified to reflect the distance used here. 
b The oxygen abundance is that derived from the total (neutral plus ion- 
ized) absorbing gas. That referring to the sole neutral gas would be some- 
how lower. 

c log(N/0) = — 1.63 ± 0.07 if the anomalously low nitrogen value in re- 
gion B4 is ignored (see Lee & Skillman 2004). 

Table 2. Average SFRs at various epochs for NGC 1705 as inferred from 
the observations (Annibali et al. 2003). 



SFR(M yr- 1 ) 


Look-back time (Myr) 


5.6 X 10~ 2 


1000-5000 


5.8 x 10~ 2 


50-1000 


7.7 x 10~ 3 


15-50 


6.8 x 10~ 2 


10-15 


0.314 


0-3 



dance in H II region B4 is ignored^ . In any case, we caution that 
nitrogen could not be determined in the same regions where direct 
electron temperatures were available. FUSE spectra seem to sug- 
gest that the N/O ratio in HI is pretty much the same, while the 
oxygen content is lower in the cold gas (Aloisi et al. 2005). 

Deep and accurate photometry has allowed to resolve stars 
from the most central regions to the extreme outskirts (Tosi et al. 
2001), so that it has been possible to divide the galaxy in several 
roughly concentric regions and to derive the SFH for each of them 
(Annibali et al. 2003): (i) different zones of the galaxy experienced 
different SFHs; (ii) the derived SF is almost continuous with fluctu- 
ations in the SFR (i.e., gasping rather than bursting); (iii) the young 
population dominates in the inner regions, where a recent, intense 
burst occurred; (iv) the strength of the most recent activity grad- 
ually decreases moving from the central to the outer regions; (v) 
NGC 1705 is definitely not a young galaxy: an old population is 
spread over the whole galaxy (see also Meurer et al. 1992); (vi) the 

t Notice that a value of log(N/0) = — 1 .63 ± 0.07 is in excellent agreement 
with the previous estimate by Storchi-Bergmann et al. (1994) - see Table 1, 
second-last row. 
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total stellar mass astrated in the galaxy is ~ 3 x 10 Mq and ~ 20 
per cent of the stars are younger than 1 Gyr; (vii) a Salpeter (a; = 
1.35) or slightly steeper (x = 1.6) field star IMF provides the best 
agreement with the observations. 

In Table 1 we summarize some structural parameters and 
observed properties of NGC 1705. Listed are: the distance to 
NGC 1705, D; its Holmberg radius, ife; its total atomic gas mass 
(HI + He), M gas ; its stellar mass, M st ar S ; its present-day metallic- 
ity, Z, and the mean oxygen, nitrogen, iron and neon abundances 
measured in H II regions and neutral gas. Table 2 reports the aver- 
age SFRs in representative age bins for the whole surveyed area, 
~ 12.5 kpc 2 , as listed in Annibali et al. (2003 - their table 6, last 
row). 

2.2.2 NGC 1569 

NGC 1569 is an exceptionally active DIG, which contains three 
SSCs (De Marchi et al. 1997; Origlia et al. 2001; Sirianni et al. 
2005) and many H II regions out of the SSC zone, where the SF is 
active now (Waller 1991). According to Aloisi et al. (2001), young 
stars (m > 8 M , r m < 50 Myr) are mostly clustered around the 
three SSCs, whereas intermediate-age objects (1.9 < m/M Q < 8, 
50 Myr < r m < 1 Gyr) are more evenly distributed. Older stars 
(m < 1.9 M Q , r m > 1 Gyr) are visible at the outskirts of the star- 
bursting region. Filaments are found at different wavelengths and 
a strong spatial correlation is present between the extended X-ray 
emission and the Ha filaments detected in the optical (Martin et 
al. 2002). In HI, a dense, clumpy ridge distribution appears, sur- 
rounded by more extended, diffuse neutral hydrogen. Discrete fea- 
tures such as arms and bridges are observed (Stil & Israel 2002) 
and the dust properties (Galliano et al. 2003) are consistent with 
the presence of shocks. Perhaps, the galaxy is presently ingest- 
ing a companion gaseous cloud (Stil & Israel 1998; Miihle et al. 
2005), and it is tempting to speculate upon the enhanced SF to co- 
incide with the accretion episode. Recent SN explosions triggered 
a GW (Martin et al. 2002), but it is unlikely that it will efficiently 
remove the ISM, unless stripping or other environmental factors 
intervene (D'Ercole & Brighenti 1999; see also Mac Low & Fer- 
rara 1999). Chandra observations reveal large inhomogeneities in 
the metal abundance of the ISM, with Z ranging from 0.1 to 1 Zq 
(Martin et al. 2002), again suggesting that multiple SN explosions 
from recent starbursts polluted the ISM. 

The SFH of NGC 1569 has been extensively studied in the last 
decade. Vallenari & Bomans (1996), from V and I optical bands, 
find that a global SF episode occurred from 100 to 4 Myr ago. They 
find hints for an older episode (from 1.5 Gyr to 150 Myr ago) and 
rule out the existence of long-lasting quiescent phases in the last 
1.5 Gyr. Similarly, based on B and V bands, Greggio et al. (1998) 
derive for the most recent burst - which, in the examined region, 
ended from 5 to 10 Myr ago - a duration of > 100 Myr and a very 
high rate, ~ 0.5 M Q yr _1 in the inner ~ 0.14 kpc 2 of the galaxy if 
Salpeter's IMF is assumed over the 0.1-120 M stellar mass range. 
If quiescent periods occurred, they lasted less than ~ 10 Myr. By 
using NIR data, Angeretti et al. (2005) looked deeper backward in 
time and found three major episodes of SF in the last 1 (possibly 
2) Gyr. They are able to distinguish two episodes at least for the 
most recent SF activity. The first started ~ 150 Myr ago and ended 
~ 40 Myr ago. The second formed stars between 37 and 13 Myr 
ago, at a mean rate about 3 times higher than the preceding one, i.e. 
~ 0.13 M Q yr _1 if Salpeter's IMF is assumed, in agreement with 
previous literature data. From their Hubble Space Telescope (HST) 
NICMOS/NIC2 data, the same authors also find evidence for a gap 



Table 3. NGC 1569 identity card. 



Quantity 


Observed value 


References 


D 


2.2 ± 0.6 Mpc a 


1 


d 


1.85 kpc 


2 


M gas 


(1.5 ± 0.3) x 10 8 M Q 


1 


Mtot 


3.3 x 10 8 M Q 


1 


h l 


0.46 ± 0.09 


1 




0.264-0.73 


3,4 


Z 


0.004 


5 


(He/H) 


0.080 ± 0.003 b 


6 


log(0/H) + 12 


8.19-8.37 


6, 7,8 


log(N/0) 


-1.39 ±0.05 


6 




-0.81 


8 


SFR 


~ 0.5 M yr- lc 


9 




0.32 M yr- 1 


10 



1- Israel 1988; 2- Stil & Israel 2002; 3- Lee et al. 2003a; 4- Martin et al. 
2002; 5- Gonzalez Delgado et al. 1997; 6- Kobulnicky & Skillman 1997; 
7- Martin 1997; 8- Storchi-Bergmann et al. 1994; 9- Greggio et al. 1998; 
10- Hunter & Elmegreen 2004. 

a Recently, Makarova & Karachentsev (2003) have derived from the I 
magnitude of the tip of the RGB two possible values for the distance to 
NGC 1569, 1.95 ± 0.2 Mpc and 2.8 ± 0.2 Mpc. 

b The error bar is likely to be larger (see Olive & Skillman 2004 reanalysis 
of helium abundance determinations in extragalactic H n regions). 
c The tabulated value refers to a look-back time of 0.1-0.15 Gyr. It has 
been derived for the central ~ 0. 14 kpc 2 of the galaxy, assuming a Salpeter 
IMF. The data also indicate a halt in the SFR 5-10 Myr ago. Strong SF 
activity is expected between 1.5 and 0.15 Gyr ago. 



in the SF activity (or, alternatively, for a low-level SF) from 300 
to 150 Myr ago, and traces of an older SF episode, occurred about 
1-2 Gyr ago at a rate between 0.01 and 0.06 Mq yr -1 (depending 
on the assumed duration). The relatively recent, strong SF activity 
proceeded at rates that could have not been sustained over a Hubble 
time, since it would have consumed all the available gas in ~ 1 Gyr, 
or even less (Greggio et al. 1998). Although a low-level SF is not 
excluded in the past, certainly the last 1-2 Gyr have been peculiar 
in the SFH of NGC 1569. There is an overall agreement among the 
various SFHs inferred from different bands. 

In a recent work, Anders et al. (2004) from multi-colour 
archive HST data show that a high level of star cluster formation 
started ~ 25 Myr ago in NGC 1569, while a lower level of activ- 
ity characterized older epochs, with a secondary peak ~ 100 Myr 
ago. The age distribution of star clusters displays a gap in the age 
interval 160^-00 Myr. Heckman et al. (1995) find that observations 
of the expanding superbubbles in NGC 1569 can be fairly well ex- 
plained assuming a SF episode at a constant rate over the ages 13- 
32 Myr, which are very similar to those inferred by Angeretti et al. 
(2005) for their youngest burst. Hunter & Elmegreen (2004) find a 
present-day SFR ip(t novj ) ~ 0.32 M Q yr -1 from Ha imaging. 

For reference, in Table 3 we summarize some basic ob- 
servational properties of NGC 1569. Listed are: the distance to 
NGC 1569, D; its maximum optical size, d; its total atomic gas 
mass (HI + He), M gas ; its total mass, M to t \ the ratio between the 
latter two, fi, which shows that from 1/4 to 3/4 of the total mass 
eventually resides in neutral atomic gas, depending on different au- 
thors; the total metal, helium, oxygen and nitrogen abundances in 
the ISM of NGC 1569 and, finally, a lower limit to the intensity of 
the latest SF episode, according to optical data. It is worth notic- 
ing that the N/O ratio measured by Kobulnicky & Skillman (1997) 
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for NGC 1569 is significantly lower than that derived by Storchi- 
Bergmann et al. (1994). We will compare only the more recent, 
more accurate determinations to our model predictions. 



3 MODELS 

In order to study the formation and evolution of late-type dwarfs, 
we use an updated version of the galactic chemical evolution (GCE) 
model developed by Bradamante, Matteucci & D'Ercole (1998). In 
particular: (i) we adopt up-to-date stellar yields, (ii) we revise the 
prescriptions about the GWs, and (iii) we treat the infall of pri- 
mordial gas as a free parameter. The adopted SFH and IMF are 
those inferred from previous CMD analyses, less prone to ad hoc 
assumptions. 

The main features of the model are the following: 

(i) the model is one zone with instantaneous and complete mix- 
ing of gas inside it; 

(ii) the instantaneous recycling approximation is relaxed, i.e. the 
stellar lifetimes are taken into account in detail; 

(iii) GWs originate when the thermal energy of the gas equates 
its binding energy; 

(iv) the winds preferably expel metals from the galaxy, while 
retaining large fractions of hydrogen; 

(v) the GW efficiency varies with time according to the SN (II 
plus la) rate predicted by the model. 



3.1 General prescriptions 

Homogeneity and instantaneous mixing of gas are reasonable 
working hypotheses, as long as gradients are not observed in the 
galaxy and the cooling time of the gas proves to be short enough. 
Hydrodynamical simulations of the behaviour of the ISM and the 
metals ejected by massive stars after a starburst suggest that most 
of the metals effectively cool off in a few million years (Hensler 
et al. 2004; Recchi et al. 2001)"^ . The existence of a small age 
gap between previous SF activity and the latest ongoing burst in 
NGC 1705 (~ 6 million years; see Table 2) is consistent with quite 
short time scales for gas cooling as well. Yet, the low occurrence 
of localized chemical pollution in the vicinity of young star clus- 
ters in star-forming dwarfs (see Sect. 2.1.1) points to time scales 
for metals cooling longer than ~ 10 7 yr (Kobulnicky & Skillman 
1997). 

Galactic outflows develop when the thermal energy of the gas, 
Eth, equates its binding energy, Eh- Following Bradamante et al. 
(1998), the thermal energy of the gas grows as a consequence of SF 
activity due to energy injection from stellar winds and SN explo- 
sions (both Type II and Type la SNe are considered): 

-Eth(i) = ^snii-Bsn / Rsmi(t') dt' + 
Jo 



3N / 
Jo 



+ »?SNIa-EsN / -RsNIa(t') di ' + 



~t~ ?7wind-E\vind 



JO J 12 



(fi(m)ip(t') dm dt' . 



(1) 



t Notice that these results strongly rest on the assumed (low) heating effi- 
ciency of SNell, i.e. on the (low) fraction of the SN explosion energy which 
remains effectively stored in the ISM and is not radiated away, jjsnii = 
0.03. The larger 7JSNII, the longer the cooling time of the ejecta. 



tp(m) is the IMF, ip(t) is the SFR, E S n = 10 51 erg is the to- 
tal energy released by a typical SN explosion, Swind = 10 49 
erg is the total energy injected into the ISM by a typical massive 
star through stellar winds. HsNii(i) and flsNia(i) are the rates 
of SN explosions, which are computed according to Matteucci & 
Greggio (1986). rysNii, »?SNia and r? W md represent the efficiencies 
with which the energies of the stellar explosions and stellar winds 
are transferred into the ISM. According to Recchi et al.'s (2001, 
2004) chemodynamical results for IZw 18, the efficiencies of en- 
ergy transfer for Type II SNe, Type la SNe and stellar winds are 
?7snii = 0.03, rysNia = 1.0 and ?7wind = 0.03, respectively (since 
SNela explode in a hotter medium, rarefied by previous SNII ex- 
plosions, their energy is more efficiently thermalized). However, in 
case of protracted SF, it is likely that SNe of both types meet with 
the same ISM conditions. Therefore, tjsnii — ?7SNia would be a 
more reasonable choice. We test both hypotheses in our study. The 
binding energy of the gas is a function of the assumed dark matter 
(DM) content and distribution. A massive (M d /Mi um = 10) dif- 
fuse (R c s/Rd = 0.1) dark halo is assumed, where Md and Mi um 
are the dark and luminous mass, respectively, and R a s and Rd are 
the effective radius and the radius of the DM core (see Table 4). 
Once the conditions for the onset of a galactic-scale outflow are 
met, metals are assumed to be lost more easily than the neutral hy- 
drogen. This is expected on the basis of the geometry and mass of 
the system: flattened galaxies of 10 8 -10 9 Mq loose preferentially 
metal-enriched matter along the polar direction, while the unpro- 
cessed gas remains almost unperturbed (De Young & Gallagher 
1990; Mac Low & Ferrara 1999; D'Ercole & Brighenti 1999; Rec- 
chi et al. 2001; Marcolini, Brighenti & D'Ercole 2004). Most of 
the metal-rich ejecta are pushed to large distances from the galaxy 
(of the order of several kpc), but their final fate is uncertain: they 
could either cool and fall back on to the galaxy plane, or be re- 
moved from the galaxy. Hydrodynamical simulations of the inter- 
action between SN-powered gas outflows and the local intergalactic 
medium (IGM) in fiducial galaxies resembling dwarf spheroidals 
show that ram pressure and tidal stripping can efficiently pull away 
matter loosely bound to the galaxy (Murakami & Babul 1999). If 
a flattened rather than a spherical ISM distribution is assumed, the 
fraction of metals retained by a moving galaxy is found to be up to 
three times larger than that retained by a galaxy at rest (Marcolini 
et al. 2004). This is interesting in the light of the empirical evidence 
that dwarf spheroidals are depleted of their gas, whereas DIGs and 
BCDs still show large gaseous contents. 

We follow in detail the time evolution of several chemical 
species (H, D, He, C, N, O, Mg, Si, S, Ca and Fe), by means of 
the following basic equation: 



( (t) dS° ut (i) 



dt 



dt 



dt 



(2) 



where Qi(t) — Xi(t)M gas (t) /M nor is the fractional gas mass in 
the form of the element i at the time t, normalized to a fixed total 
mass, M nor . The quantity -X»(t) represents the abundance by mass 
of the element i at the time t and the summation over all the ele- 
ments in the gas mixture is equal to unity. ip(t) is the SFR, while 
1Zi(t) represents the fraction of mass which is restored into the ISM 
in the form of the element i by dying stars. A detailed description of 
this term, which contains all the prescriptions about the stellar IMF, 
stellar nucleosynthesis and Type la SN progenitors, can be found in 
Matteucci & Greggio (1986). The specific yield sets adopted in this 
work are discussed in Sect. 3.2. Finally, the last two terms on the 
right-hand side of Equation (2) account for possible gas inflows and 
outflows, respectively. 
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Table 4. Parameters of the models with their adopted values. 



Parameter 


Meaning 


Adopted value 
NGC 1705 NGC 1569 




normalizing mass (gaseous matter available for accretion inside the DM 
potential) 


7 x 10 8 M 5 x 10 8 M Q 


M d 


DM amount 


7 x 10 9 M 5 x 10 9 M 


-Rcff 


effective radius of the visible matter 


1.2 kpc 1 kpc 




effective to DM core radii ratio 


0.1 


T 


time scale for mass accretion 


8 Gyr 0.5, 5 Gyr 


^SNII 


SNII thermalization efficiency 


0.20 0.03, 0.20, 0.50 


»?SNIa 


SNIa thermalization efficiency 


0.20 0.20,0.50,1.00 


^?wind 


stellar wind thermalization efficiency 


0.20 


tnow 


age of the universe 


13.7 Gyr 



In many literature models for dwarf galaxy formation, DIGs 
and BCDs form from continuous infall of gas of primordial chemi- 
cal composition, until a mass from a few 10 8 M up to ~ 10 9 M 
is accumulated. The rate of gas infall is normally parametrized as 



dt 



Al nf e- t/T 



(3) 



T (l_ e -tnow/r)' 

with r, the infall time scale, being the most critical free parameter 
and largely varying amongst different authors (e.g., r = 0.5 Gyr for 
Matteucci & Tosi's 1985 models; r = 1, 4, 7 Gyr for Mouhcine & 
Contini's 2002 models). Here we run models with the choices for 
r specified in Table 4. X\ ni — Xf is the primordial abundance of 
element i in the inflowing gas. Here we adopt (D/H)p = 2.5 x 10~ 5 , 
( 3 He/H) P = 0.9 x 10~ 5 , and Y P = 0.248 (see Romano et al. 2003, 
and references therein). The quantity i n ow, the present time, is set 
to 13.7 Gyr. In addition, we test also the cases of a constant infall 
rate, 



= const, 



dg| nf (t) 
dt 

and of an increasing exponential, 

dGr ! ( t) _ X! nf e* /T 
dt 



(4) 



(5) 



r ( e i„ow/r _ iy 

with the lowest accretion rates at the early stages of galaxy for- 
mation and the highest ones at later epochs, to coincide with 
the observed late SF boosting. Notice that the infall rate is 
normalized to the quantity M nor , which obeys the condition: 

IK AV IK >YV now )■, 

where M gaa (t now ), M stars (i now ), and M rom (t now ) are the masses 
in form of gas, stars and stellar remnants, respectively, and 
A^Gw(inow) is the total mass ejected from the galaxy by GWs at 
the present time. 

The rate of gas loss via GWs at the time t is assumed to be 
proportional to the amount of gas present at that time. For each 
element i, 



dt 



■■ Wi {t)xr\t)Q(t), 



(6) 



where X° ut (t) = Xi (t) . Here we call for a dependence of Wi (t) 
on time via a dependence on the Type II plus Type la SN rate (ex- 
pressed in number Gyr -1 ): 



Wi(t) = /3i[r7SNII-RsNIl(i) + 7?SNIa-RsNIa(t)]. 



(7) 



The parameter $ is a proportionality constant indicating the effi- 
ciency of mass loss for a specific element i. In the case of differen- 
tial winds, the value of j3% is assumed to be different for different 



elements. In particular, higher f)i values are assigned to heavy ele- 
ments produced in SNe, while only a minor fraction of the neutral 
hydrogen is allowed to leave the galaxy. 
The SFR has the simple form: 



V(t) = vQ(t), 



(8) 



where v, the efficiency of SF, is expressed in units of Gyr -1 , and is 
fixed so as to recover the SFR inferred from the observations. 

The adopted IMF (by mass) is an extension of that originally 
proposed by Salpeter (1955) to the whole stellar mass range: 



<p(m) = Cm 



(9) 



with x = 1.35 for 0.1 < m/M < 100. C is the normalization 
constant, and is obtained by imposing the condition 

too 

(p(m)dm = 1. (10) 



I 



This IMF is in agreement with the results of the CMD analysis for 
both NGC 1705 and NGC 1569, though one must be aware that the 
synthetic CMDs actually constrain the IMF slope of these galaxies 
only for m 6 M (Annibali et al. 2003). Flatter slopes in the low- 
mass range cannot be excluded (Angeretti et al. 2005). We discuss 
the effect on model results of adopting different IMF formulations 
- still consistent with the CMD analysis - in Sect. 4.3.1. 

3.2 Nucleosynthesis prescriptions 

To quantitatively estimate the uncertainty on the results related to 
different nucleosynthesis assumptions, we have run our models as- 
suming various sets of stellar yields. 

For low- and intermediate-mass stars (LIMS), i.e., stars with 
0.8 < m/M < 8, both the yields of van den Hoek & Groenewe- 
gen (1997) and those computed by Meynet & Maeder (2002) for 
rotating stellar models are considered. For the van den Hoek & 
Groenewegen (1997) yields, we adopt either those computed with 
constant mass loss parameter along the asymptotic giant branch 
(AGB), 7/agb = 4 (Table 5; models labelled N), or those computed 
with metallicity-dependent ijagb (Table 5; models labelled W), or 
those computed with the minimum hot bottom burning (HBB) ex- 
tent allowed by the observations (Table 5; models labelled H). For 
the Meynet & Maeder (2002) yields, we adopt those computed for 
initial rotational velocities of the stars of Vi n { = 300 km s _1 (Ta- 
ble 5; models labelled M). 

For massive stars, i.e., stars with m > 8 M , the adopted 
nucleosynthesis prescriptions are either from Nomoto et al. (1997; 
Table 5; models N and H), or from Woosley & Weaver (1995; 
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Table 5. SFH, GW efficiency and nucleosynthesis prescriptions. 



Star formation history 

Model Number of bursts Time of burst occurrence Duration of the burst v 



(Gyr) (Gyr) (Gyr' 1 ) 

1 Star formation history from Annibali et al. (2003) 

2 3 5.55/11.55/13.6 6/2/0.1 0.05/0.15/2 

3 3 5.55/12.55/13.6 6/1/0.1 0.05/0.3/3 

4 5 4.575/6.575/8.575/11.55/13.6 0.25/0.25/0.25/2/0.1 0.5/0.5/0.5/0.15/2 

5 3 0.7/11.55/13.6 10.5/2/0.1 0.012/0.15/2.4 



GW ejection efficiencies (ft) a 

H, D, He C, N, O, and heavier 

a 2 X Kr b 2 X 10- b 

b 2 x Kr 6 2.5 x 10~ 5 
c As above, except for nitrogen, for which /3n = 5.5 X 10 -5 

d 0.5 X 10 _s 2.5 X 10 -5 

e 0.5 x 10" 5 1.5 x 10" 5 

f As case d, except for nitrogen, for which = 5.5 X 10 — 5 



Nucleosynthesis prescriptions 
LIMS Massive stars 

N van den Hoek & Groenewegen 1997, ijagb const Nomoto et al. 1997 D 

W van den Hoek & Groenewegen 1997, »7agb var Woosley & Weaver 1995 b 

H van den Hoek & Groenewegen 1997, minimum HBB Nomoto et al. 1997 b 

M Meynet & Maeder 2002 c Meynet & Maeder 2002 c 



a Notice that the actual GW ejection efficiencies are given by Wi(t) = ft[»7sNll-RsNll(*) + »?SNla-RsNla(t)], where -RsniiW and ^?SNla(*) are the 
Type II and Type la SN rates by number at time t, respectively, each one multiplied by the appropriate SN heating efficiency. 

b The yields of 12 C for stars more massive than 40 Mq were multiplied by a factor of 3 in order to fit the Milky Way carbon abundance data (Chiappini et 
al. 2003a). 

c We adopt the set of yields computed for initial rotational velocities of the stars of t)j n j = 300 km s -1 . 



Table 5; models labelled W), or from Meynet & Maeder (2002; 
Table 5; models labelled M). 

Type la SN nucleosynthesis is accounted for by adopting the 
prescriptions by Thielemann, Nomoto & Hashimoto (1993), their 
model W7. 

Each model is characterized by one number and one pair of 
letters; the number refers to the adopted SFH, the small letter to the 
adopted GW efficiency (which, according to the differential wind 
picture, is higher for the heavy elements than for the light ones), 
the capital letter to the adopted stellar yields (see Table 5). 



4 RESULTS 

Most GCE studies published so far assume that DIGs and BCDs 
evolve through successive bursts of SF - which possibly induce 
GWs - alternating long quiescency periods (Matteucci & Chiosi 
1983; Matteucci & Tosi 1985; Pilyugin 1993; Marconi et al. 1994; 
Bradamante et al. 1998). Very few attempts are made at modelling 
specific objects (e.g., Carigi et al. 1999; Recchi et al. 2001, 2004, 
2005). However, given the high number of free parameters involved 
in dwarf galaxy modelling, it is instead of fundamental importance 
to model individual galaxies whose SFH and IMF have been de- 
rived from observations. 



4.1 NGC 1705 

NGC 1705 has its SFH and IMF slope known thanks to deep HST 
photometry probing epochs up to 5 Gyr ago (see Sect. 2.2.1). The 
stellar mass astrated at the inferred SFR over the last 5 Gyrs is suf- 
ficient to account for the stellar mass we observe in the galaxy to- 



day (Annibali et al. 2003). Therefore, there is no need to speculate 
upon earlier SF activity, and we can safely assume that the galaxy 
formed the vast majority of its stars during the last ~ 5 Gyrs at the 
average rates suggested by Annibali et al. (2003; see Table 2). Fol- 
lowing those authors, the stellar masses are distributed according 
to a Salpeter IMF (x = 1.35) over the whole (0.1-100 Mq) stellar 
mass range, unless otherwise stated. 

In Fig. 1 we show the adopted SFH together with some results 
from Model iaiV. The adopted structural parameters of the model 
are: M nor =7 X 10 8 Mq, R cti = 1.2 kpc, M d = 7 X 10 9 Mq, 
Reff/Rd = 0.1 (see Sect. 3.1 and Table 4 for explanations regard- 
ing model parameters). The assumed infall law is exponentially de- 
creasing with time with e-folding time r = 8 Gyr (Fig. If, dashed 
line). This is a quite long time scale, comparable to that for the for- 
mation of the solar neighbourhood (tsn = 7 Gyr; Chiappini, Mat- 
teucci & Gratton 1997). It allows a non-negligible fraction of the 
cold gas to be available for SF at late times. The assumed thermal- 
ization efficiencies from SNe and stellar winds are 7)snii = »?SNia = 
Vwind = 0.20. 

Model laN well reproduces the present-day gaseous (Fig. la) 
and stellar (Fig. lc) contents of NGC 1705, as well as its current 
metallicity (Fig. Id). In the framework of this model, a galactic 
outflow develops very soon, already ~ 250 million years after the 
onset of the SF, due to the strong SN feedback. We predict a current 
mass loss rate of ~ 0.04 Mq yr _1 through the outflow, i.e., about 
an order of magnitude lower than the present-day SFR. This is due 
to the fact that SNe with progenitors born during the last sudden 
burst of SF, have not had time to explode and leave their signature 
on the GW properties yet. Therefore, it is still the previous SF ac- 
tivity which determines the features of the outflow in our model. An 
attempt to evaluate the actual mass loss rate of NGC 1705 has been 
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Figure 1. a) M gas , b) SFR, c) M s t a rs, d) gas metallicity, e) Type II (solid line) and Type la (dashed line; multiplied by a factor of 10 to make it clearly visible) 
SN rate, f) M GW (solid line) and M inf (dashed line) as functions of time; g) [N/Fe] and h) [0,Fe] as functions of [Fe/H] for Model laN for NGC 1705. The 
stars represent observed values. Both the theoretical and the observed (Aloisi et al. 2005) [el/Fe], [Fe/H] ratios are normalized to solar values by Asplund et 
al. (2004). 



made by Meurer et al. (1992) in their multi-wavelength study of this 
object. They parameterize the rate of mass loss from the galaxy in 
terms of two quantities, Ma ow , the total mass of gas entrained in 
the outflow, and M to t /Lb, the total mass-to-light ratio. The uncer- 
tainty in Mfl ow (nearly three orders of magnitude) and the poorly 
known DM contribution to M to t make the observed mass loss rate 
really ill-constrained: Meurer et al. indicate a minimum mass loss 
rate of 0.0026 M© yr" 1 [one order of magnitude below our predic- 
tion, MGw(tnow) — 0.04 Mq yr -1 ], but values up to three orders 
of magnitude higher are also permitted. It is thus clear that theo- 
retical models for the development of GWs are poorly constrained 
by the observations currently available for NGC 1705. In this re- 
spect, a much more interesting case to study is NGC 1569, which 
has its outflow much better characterized by deep Chandra spec- 
tral imaging (Martin et al. 2002). Therefore, we defer any further 
consideration on the GW behaviour to Sect. 4.2, where we treat in 
detail the formation and evolution of NGC 1569. 



As mentioned above, we need to assume that the winds driving 
the evolution of NGC 1705 are differential, i.e., they carry away a 
large fraction of the newly produced metals, while retaining most of 
hydrogen. This is necessary in order to reproduce both the observed 
present-day gas content and the current metallicity of the galaxy. In 
Model laN, we set the f3i parameters to 2 x 10 -5 for the heavy 
elements and to 2 x 10 -6 for the light ones (see Sect. 3.1 and Ta- 
ble 5). Setting f3i to 2 x 10~ 6 for all the chemical species would 
result in a present-day metallicity of 0.017, i.e. almost solar. On 
the contrary, setting /3; to 2 x 10~ 5 for all the chemical elements 
would result in a present-day metallicity of 0.04, i.e. twice solar, 
since in this case the metals are mixed in a much lower gaseous 
mass: Mg as (t now ) ~3 x 10 6 Mq, at variance with the observa- 
tions. Clearly, the model is not unique, and the formulation for the 
Wi(t) parameters is allowed to change. However, in agreement with 
previous analysis (Pilyugin 1993; Marconi et al. 1994; Recchi et al. 
2001, 2004, 2005) we find that models assuming a common value 
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for Wi(t) (i.e., no differential winds) unavoidably fail to reproduce 
(at least some of) the observational constraints. 

4.1.1 Detailed chemical properties of the neutral and ionized 
gas phases in NGC 1705 

The abundances of oxygen and nitrogen in the local ISM of 
NGC 1705 have been derived from spectroscopy of HII regions 
by several groups (see Sect. 2.2.1). In the following, we adopt 
for comparison with our model predictions the up-to-date val- 
ues of log(0/H)+12 and log(N/0) suggested by Lee & Skillman 
(2004 - Fig. 2, filled diamonds, with 1- and 2-cr error bars). While 
the present-day oxygen abundance is very well reproduced by the 
model (Fig. 2, upper panel, cfr. the final points of the tracks), the 
current N/O ratio is largely overestimated (Fig. 2, lower panel). In 
Fig. 2 we show results obtained by adopting either x = 1.35 over 
the whole stellar mass range (solid lines) or x = 1.6 above 2 Mq 
(dashed lines; Annibali et al. 2003). From now on, we will always 
refer to models for NGC 1705 as models computed with x = 1.35 
over the whole stellar mass range. 

Recently, FUSE and Space Telescope Imaging Spectrograph 
(STIS) Echelle spectra of NGC 1705 have allowed the determina- 
tion of the metal (O, Ar, Si, Mg, Al, N, and Fe) abundances also 
in the neutral medium of NGC 1705 (Aloisi et al. 2005). However, 
when comparing them to model predictions, one must be aware 
that saturation (especially for O I), photoionization, contamination 
by gas laying along the line of sight, depletion into dust grains (e.g., 
Fe could be more easily locked into grains compared to O) all could 
affect the derivation of the abundances from the neutral gas (Aloisi 
et al. 2005). In Fig. lg,h we display [N/Fe] and [O/Fe] as functions 
of [Fe/H]. All the ratios are normalized to the solar values of As- 
plund, Grevesse & Sauval (2004). At [Fe/H] ~ -1.2, the [N/Fe] 
predicted by Model JaN is higher than the one measured in the 
neutral gas. This is not unexpected on the basis of the results for 
the ionized gas phase shown in Fig. 2. The [O/Fe] ratio is overes- 
timated by the model as well. It is worth emphasizing that, while 
the observed [N/H] and [Fe/H] were corrected for the presence of 
ionized gas lying along the line of sight, [O/H] was not (Aloisi et 
al. 2005). 

The pronounced discrepancy between model predictions and 
observations for nitrogen - the predicted nitrogen abundance is 
sistematically higher than observed, both for the neutral and the 
ionized gas phases - requires further analysis, especially since the 
model well reproduces the other reliable observational constraints. 
In the remaining part of this section, we discuss how current uncer- 
tainties in the theories of stellar evolution and nucleosynthesis may 
affect our results. 

4.1.1.1 Less hot bottom burning in intermediate-mass stars? 

In Fig. 3 we display the evolution of [N/Fe], [O/Fe] vs. [Fe/H] and 
log(0/H)+12, log(N/0) vs. time predicted by models for NGC 1705 
adopting one of the following combinations: (i) the yields from 
massive stars by Nomoto et al. (1997) plus the yields from LIMS 
by van den Hoek & Groenewegen (1997), either with a mass loss 
parameter along the AGB independent of metallicity - Models JaN 
(first column, solid line) and lcN (first column, dashed line) - or 
with the minimum HBB extent (Model JaH, second column); (ii) 
the yields from massive stars by Woosley & Weaver (1995) plus the 
yields from LIMS by van den Hoek & Groenewegen (1997) com- 
puted with metallicity-dependent mass loss parameter along the 
AGB (Model laW, third column); (iii) stellar yields by Meynet & 
Maeder (2002) for rotating stellar models (Models JaM and JcM, 
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Figure 2. Temporal evolution of oxygen f top panel ) and nitrogen to oxy- 
gen (bottom panel) in the ISM of NGC 1705. Predictions from Model JaN 
(solid lines: Salpeter IMF; dashed lines: steeper IMF - see text for details) 
are compared to the available observations. The present-day oxygen abun- 
dance of NGC 1705 is that given by Lee & Skillman (2004; mean value of 
[O m] A 4363 measurements; small filled diamond). The oxygen abundance 
of the neutral ISM is that derived from far-UV observations (Aloisi et al. 
2005; dashed horizontal hand), not corrected for contamination by ionized 
gas lying along the line of sight. The current N/O ratio in NGC 1705 is that 
suggested by Lee & Skillman (2004; small filled diamonds), with (lower 
value) and without (higher value) region B4. A value of log(N/0) = —1.82 
± 0.07 characterizes the neutral phase (Aloisi et al. 2005; dashed horizontal 
hand), which is the same as HII regions within the errors. 



fourth column, solid and dashed lines, respectively). We recall that 
Meynet & Maeder (2002) have not computed the third dredge-up 
and HBB phases; hence their stellar models for LIMS are somewhat 
incomplete. Model JcN (lcM) differs from Model JaN (JaM) 
only in the adopted values of the GW efficiencies. Models JcN and 
JcM are discussed in detail in Sect. 4.1.2. 

Uncertainties in the physics of stellar models for intermediate- 
mass stars cause the predicted log(N/0) at the present time to vary 
widely, from -0.9 (Model JaW) to -1.5 (Model JaM). Stellar 
models with rotation and without HBB give rise to the lowest theo- 
retical N/O ratios, but they are still only marginally compatible with 
the observations. Results on the oxygen behaviour are firmer. Yet, 
lower [O/Fe] ratios are predicted during the whole galaxy evolu- 
tion when adopting the yields of Woosley & Weaver (1995) rather 
than those by Nomoto et al. (1997) for massive stars. Meynet & 
Maeder's yields are the only ones able to reproduce both the ob- 
served O/H and N/O at the present time. However, the discrepancies 
among model predictions and observations for the neutral phase, 
though sensibly reduced, are not eliminated (Fig. 3, last column, 
solid lines). 

4.1.1.2 Localized massive star pollution in NGC 1705? 

Notwithstanding the large uncertainties affecting the nitrogen 
yields, the discrepancy between model predictions about nitrogen 
evolution in NGC 1705 and the relevant observations is hardly ex- 
plained only in terms of the adopted stellar yields. A second pos- 
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Figure 3. [N/Fe] and [O/Fe] as a function of [Fe/H], log(0/H)+12 and log(N/0) as a function of time. Predictions from Models laN, laH, laW and laM 
(solid lines) as well as Models IcN and lcM (dashed lines) are displayed for comparison one to another and with data for both the neutral (Aloisi et al. 2005; 
stars and dashed strips) and ionized gas (Lee & Skillman 2004; small filled diamonds with 1- and 2-cr error bars). 



sibility is localized self-pollution from dying young massive stars 
born during the last 3 Myr of SF activity. According to Meynet & 
Maeder (2002), the N/O ratio in the ejecta of a 60 Mq star, having a 
lifetime of ~ 3.5 million years, is about log(N/0) ~ —2.2 for Zi n i = 
0.004. It would be more interesting to compare to the observations 
the log(N/0) values predicted for rotating massive stars with life- 
times shorter than 3 million years, but, unfortunately, no star mod- 
els with initial mass larger than 60 Mq are provided yet. Anyway, 
we notice that the N/O ratio in the ejecta of massive stars with ini- 
tial metallicity Zi n i = 0.004 and initial rotational velocity Vi n j = 
300 km s _1 is a decreasing function of mass (Meynet & Maeder 
2002). Therefore, lower N/O ratios are expected in case of localized 
massive star self-pollution with respect to the general predictions 
of homogeneous GCE models, which may solve the discrepancy 
between model predictions and observations. However, such a so- 
lution would not account for the low nitrogen abundance measured 
with FUSE. Moreover, the constancy of log(0/H)+12 among dif- 
ferent H II regions issues a serious challenge to this interpretation. 



4.1.2 The role of gas flows 

IMF variations in the range allowed by CMD analyses do not al- 
low to solve the 'nitrogen problem' in NGC 1705 (see Fig. 2, lower 
panel). We are thus left with the possibility that the efficiency of 
nitrogen ejection through the wind of NGC 1705 is higher than 
that for oxygen, i.e., /3n//3o 1. Though this hypothesis needs 
to be tested by detailed hydrodynamical computations before it can 
be accepted or rejected, we compute two models for NGC 1705 
in which a higher /3n//?o value is assumed, namely, /3n//3o — 
2 rather than 1 (Models IcN and JcM), and find that the latter 
matches fairly well the available observations (Fig. 3, last column, 
dashed lines). In conclusion, it seems that the outflowing gas should 
be differentially enriched by the products of stellar nucleosynthesis 
in order to make the model predictions match the observations. 

At this point, one must check if our conclusions change by 
modifying the prescriptions on gas inflows. Up to now, we assumed 
that some unprocessed gas is continuosly accreted by the galaxy 
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Figure 4. Same as Fig. 1, but for models for NGC 1705 where the infall rate is either constant (thick solid lines) or exponentially increasing in time (thick and 
thin dashed lines). In the latter case, the hypothesis of differential GWs may be relaxed (see text for details). 



during the whole Hubble time according to an infall law exponen- 
tially decreasing with time§ and with a time scale r = 8 Gyr. In the 
following, we discuss the results obtained when considering an in- 
fall rate either constant or exponentially increasing with time (again 
with r = 8 Gyr). The latter case would simulate important mergers 
with gaseous lumps at late times. 

Results of models in which either a constant or an increas- 
ing infall rate are considered are sketched in Fig. 4. Once a con- 
stant infall rate is assumed, predictions from Model laN (Fig. 4, 
solid lines) almost do not differ from those obtained by the same 
model with a mass accretion rate exponentially decreasing in time 
(cfr. Fig. 1), except for the current mass loss rate, MGw(inow) = 
0.1 M yr _1 , rather than 0.04 M© yr _1 . In the case of con- 
stant infall, metals are mixed in a slightly higher gaseous mass 
[MgasOnow) = 2.1 x 10 8 M Q rather than 1.9 x 10 8 M Q ], and 

§ This is a widely adopted, but completely arbitrary, assumption of GCE 
models (see, e.g., Mouhcine & Contini 2002 and discussion therein). 



one gets Z = 0.0032 rather than Z = 0.0037 at the present time, 
a negligible difference. Instead, if the rate of mass accretion in- 
creases in time (see Fig. 4f, dashed line), the current gaseous mass 
of the galaxy is overestimated by the model, unless one does not 
raise the efficiency of hydrogen loss through the wind. In Fig. 4a 
we display the behaviour of the gas mass with time predicted by 
Models ldN and ieiV, both assuming (3i = 0.5 x 10~ 5 for the light 
elements (dashed lines). This is the highest f3i value for hydrogen 
ejection we can impose while still retaining enough mass to match 
the observations. The current metallicity of the ISM is a further 
constraint to be fulfilled, and it restricts the possible choices of f3i 
for the heavy elements. As an example, in Fig. 4d we show re- 
sults obtained with /?; equal to either 2.5 x 10 -5 (Model ldN, 
thin dashed line) or 1.5 x 10~ 5 (Model leN, thick dashed line) 
for elements heavier than helium. It is immediately seen that an in- 
fall rate exponentially increasing in time allows one to relaxe the 
differential wind hypothesis: /3hcavy/Aight = 10 for Model iaiV; 
/3heavy/Aight = 5 for Model ldN (see also Table 5). At this point, 



© 2005 RAS, MNRAS 000, 1-20 



12 D. Romano et al. 



L, 0.4 
0.3 
^ 0.0 



0.005 



0.000 



8.5 
8.0 



(A! 

' — I 

+ 

£ 7,5 

? 7.0 



"J 
o 



o 



-0.4 

-0.6 

-0.8 

-1.0 
-0.6 
-0.8 
-1.0 
-1.2 
-1.4 



2aN 






3aN 






4aN 



_M 




4aff 




5aW 





10 10 10 10 10 

t [Gyr] t [Gyr] t [Gyr] t [Gyr] t [Gyr] 



Figure 5. SFR, global metallicity and CNO abundances as functions of time according to Models 2aN, 3aN, 4aN, 4aW and 5aN for NGC 1569. See text for 
details. 



one might wonder whether a fine tuning of model parameters would 
make the Aeavy/Alght value further approach unity, eventually 
getting rid of the need for metal-enriched GWs. We ran several 
models, but we never reproduced simultaneously both the current 
gaseous mass and ISM metallicity with Aeavy /Aight = 1> unless 
making some very ad hoc assumptions with possible huge episodes 
of mass accretion at late times. However, NGC 1705's outer optical 
isophotes appear round and regular, suggesting a relatively peace- 
ful environment (Meurer et al. 1998). Finally, notice that chang- 
ing the infall law does not solve the above-mentioned discrepancy 
between model predictions and observations for nitrogen: models 
with /3n//9o = 1 never match the data, independently of the adopted 
infall rate formulation (cfr. Fig. 4g,h and Fig. lg,h). The current 
rate of mass loss predicted by models with rate of gas accretion 
exponentially increasing with time is MGw(^now) = 0.2 Mq yr _1 , 
owing to the higher hydrogen mass entrained in the outflow in this 
case. 



4.2 NGC 1569 

NGC 1569 is another well-studied, relatively nearby late-type 
dwarf, whose spectacular GW has been powered by an unusually 
high-level SF activity (see Sect. 2.2.2). Given the similarities be- 
tween NGC 1705 and NGC 1569, we find it natural to start a de- 
tailed study of NGC 1569 and see whether models adopting the 
formalism developed for NGC 1705 satisfy the observational con- 
straints available for NGC 1569 as well. 

Our models for NGC 1569 all assume that two major bursts of 
SF have occurred in the last ~ 1-2 Gyr, in agreement with the in- 
ferences from the observed CMDs (Greggio et al. 1998; Angeretti 
et al. 2005). Because of observational uncertainties and since the 
observed CMDs are relevant to small inner galactic regions - while 
we need a global SFR referring to the whole galactic area - we still 
have some degrees of freedom when choosing the actual SFRs at 
late times. Moreover, we must allow for some SF activity at ear- 
lier times, in order to build up a total stellar mass which matches 
that inferred from the observations. Possible choices are not com- 
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Figure 6. Same as Fig. 5, but for Models 5aH, 5aM. 



pletely arbitrary, since they must obey two more basic, stringent 
constraints: 

(i) Most of the latest SF activity is concentrated within the area 
surveyed by Greggio et al. (1998). Therefore, it is safe to adopt the 



SFR suggested by those authors, ~ 0.5 M Q yr , as representa- 
tive of the overall SF activity across the whole body of the galaxy. 
Besides this value, we also test the slightly lower one derived 
by Hunter & Elmegreen (2004) from Ha imaging of NGC 1569, 
~0.3M© vr" 1 . 

(ii) SF activity prior to that surveyed by HST likely proceeded 
at much lower rates (see the discussion in Angeretti et al. 2005). 
Therefore, for ages between ~ 1-2 and ~ 13 Gyr ago, we assume 
a continuous SF activity, which proceeds at rates much smaller than 
assumed for the recent episodes. However, only for the sake of 
comparison with old 'bursting' models, we also examine the case 
of a bursting SF regime at early times, with three bursts alternating 
long quiescency periods. 

Table 5 summarizes the SFHs assumed in this paper for NGC 1569 
in the light of the above considerations (Models from 2 to 5). 
Model results are discussed below. 



4.2.1 Bursting, gasping or continuous SF regimes? 

Figs. 5, 6 and Table 6 summarize and compare the results obtained 
with different models for NGC 1569 where the SFH and stellar 
yields are allowed to change one at a time. The adopted structural 
parameters of the models are: A/ nor = 5 x 10 8 M©, R c g = 1 kpc, 
M d = 5 x 10 9 M , i? cft /Rd = 0. 1 (see Sect. 3. 1 and Table 4). The 
infall law is exponentially decreasing with time with r = 0.5 Gyr 
for all the models (see, however, Sect. 4.2.3). Models 2aN and 3aN 
assume that a long-lasting SF activity preceded the strong last 1-2 
Gyr activity detected with HST. Following Greggio et al. (1998) 
and Angeretti et al. (2005), we assume that the most recent SF 
activity proceeded through two violent successive bursts: the first 
started almost 2 (Model 2aN) to 1 (Model 3aN) Gyr ago and lasted 
nearly 2 (Model 2aN) to 1 (Model 3aN) Gyr, the second started 
100 Myr ago and finished a few Myr ago. A 50 Myr gap separates 
the two episodes; assuming a longer interburst period (150 Myr - 
Angeretti et al. 2005) does not change the results. A weaker (y = 
0.05 Gyr -1 ), but long-lasting (At = 6 Gyr) SF episode is assumed 
to have occurred from t = 5.55 to t = 1 1.55 Gyr, where no detailed 
information is available from HST photometry. This pristine SF 
raises the metal content of the gas from Z = 0toZ~2x 10 -3 . For 
Models 4aN and 4a W, the most ancient activity consists of three 
intense short-lasting bursts, while Models 5aN, 5aH, 5aM adopt a 
continuous, very low-level SFR over almost the whole Hubble time. 
All the models share the same stellar nucleosynthesis prescriptions, 
but Models 4a W, 5aH, 5aM. The adopted GW efficiencies are al- 
ways the same (see Table 5). 

First, we want to emphasize that models assuming either a 
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Table 7. Model results for NGC 1569 (at t = t now = 13.7 Gyr; fourth to twelfth cols.). Model prescriptions are the same of Model 5aN (first row), except for 
the adopted SNII and SNIa thermalization efficiencies (first and second cols.) and, eventually, DM content (third col; see text for more details). The time of 
the onset of the GW, *gw. i s a l so indicated (last col.). 
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Figure 7. Distributions of low-mass, long-lived stars as functions of time 
(left panel) and metallicity (right panel) according to Models 2aN and 4aN 
for NGC 1569. The models assume the same SFH in the last ~ 2 Gyr, where 
more stringent constraints are available from HST photometry. As far as 
earlier activity is concerned, Model 2aN adopts a continuous, low-level SF 
over ~ 6 Gyr, while Model 4aN assumes three short-lasting bursts. On the 
right-hand panel the two theoretical distributions turn out to be indistin- 
guishable. 

bursting, or a gasping, or a fairly continuous, highly variable SF 
regime are all able to produce results consistent with the obser- 
vations (e.g., Models 2aN, 4aN, 5aN, Fig. 5, first, third and last 
col., respectively). They provide relatively small differences in the 
predicted frequency distributions of low-mass stars with metallic- 
ity (Fig. 7). Only much deeper CMDs could unveil the 'true' SFH 
of NGC 1569 at epochs prior to 1-2 Gyr ago. So far, we have some 
clues that early SF in NGC 1569 went on at rates significantly lower 
than the ones characterizing later episodes (Angeretti et al. 2005). 
This information, coupled to current observational evidence point- 
ing to fairly continuous SFHs for most DIGs and BCDs (Schulte- 
Ladbeck et al. 2001; Annibali et al. 2003; Tosi 2003), leads us to 
indicate the SFH labelled 5 as the most likely evolutionary scenario 
for NGC 1569. 

Second, we find that for NGC 1569 different combinations of 
stellar yields reproduce the N/O data (at 2-a, but fine tuning of 
model parameters may reconcile model predictions to observations 



even better), except for Model 4a W, adopting the van den Hoek & 
Gronewegen (1997) yields with metal-dependent t^agb for LIMS. 
This model indeed largely overestimates the N/O value measured 
in the local ISM of NGC 1569 [log(N/0) = -1.39 ± 0.05; Kobul- 
nicky & Skillman 1997]. Successful models for NGC 1569 all have 
Pn / Po = 1, at variance with models for NGC 1705, where, in order 
to fit the observations, a value of Pn/po greater than 1 is always 
required, independently of the choice of the stellar yields. 

The issue of the GW efficiencies deserves further considera- 
tions. The SF activity in NGC 1569 triggers and sustains an outflow 
on a galactic scale in our model galaxies. The more efficient is the 
SF process, the more effective is the wind in removing the newly 
produced metals, thanks to its dependence on the SN rate. For in- 
stance, Model 3aN (Fig. 5, second col), computed assuming a SF 
efficiency higher than Model 2aN (Fig. 5, first col.) during the last 
1 Gyr, does not overproduce the present-day oxygen content and 
overall metallicity of the gas even if the Pi values are kept the same, 
because the efficiency of gas removal is increased in lockstep with 
the predicted higher Type II plus Type la SN rate. In particular, we 
find that the assumed time modulation of the mass loss rate [Equa- 
tion (7)] leads to a predicted mass ejection rate at the present time 
which is of the same order of magnitude of the present-day SFR, 
in agreement with recent Chandra observations (Martin et al. 2002; 
cfr. sixth vs twelfth col. in Table 6). This would not be the case if the 
mass loss rate were simply proportional to the gas content. In that 
case, the rate of mass loss through the wind would be maximum 
when the thermal energy equals the binding energy of the gas for 
the first time, after which it would monotonically decline in time, 
owing to the lower and lower gas amount left over by the ongo- 
ing SF process. In other words, a simple dependence of the outflow 
rate on the gas mass would not allow to predict enhanced mass loss 
rates to coincide with enhanced SF episodes, as observed. It is also 
worth emphasizing that by assuming a GW simply proportional to 
the SFR (as often found in the literature) one forgets the effect of 
delayed SNIa explosions and the wind goes to zero as soon as the 
SF stops. 



4.2.2 The energetics of the ISM 

The energetics of the ISM is a complicated issue, and no recipes 
can be found in the literature which treat it without introducing a 
number of uncertain free parameters. The energy injected into the 
ISM by a typical massive star through stellar winds during its life- 
time is iJ w i n d — 10 49 erg, to be compared with an initial blast wave 
energy of £sn — 10 51 erg for both SN types. What really matters is 
the fraction of this energy which is not radiated away and remains 
effectively stored in the ISM. In previous work on IZw 18, Recchi 
et al. (2001, 2004) assigned a low thermalization efficiency of only 
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3 per cent to SNell, whereas a much higher value, 100 per cent, was 
assumed for SNela, exploding in a medium already heated by pre- 
vious SNII activity after a non-negligible time lag from the onset of 
the SF. However, those values are likely to apply strictly only to the 
case of single, instantaneous bursts of SF. If the SF is an almost con- 
tinuous process, tjsnii must probably be increased (e.g. Strickland 
& Stevens 2000). Moreover, 7/snii should vary with time, since it 
depends on the amount of gas in form of cold clouds at each time 
(Melioli & de Gouveia Dal Pino 2004). Since we do not follow the 
exchange of matter between cold and hot gas phases in our mod- 
els, we cannot state how much gas is in the form of cold clouds 
at each time. Therefore, we do not vary 7/snii and r/sNia with time, 
and assume average constant values for these parameters during the 
whole galaxy evolution. Moreover, we do not treat SF in small (pc- 
sized) zones, but we assume that the SF process is distributed over 
the whole physical dimensions of the object, i.e. on kpc scale. Yet, 
our simplistic models tell us something about the mechanisms reg- 
ulating dwarf galaxy formation and evolution. In the previous case 
of NGC 1705, we have assigned the same thermalization efficiency 
of 20 per cent to stellar winds, SNell and SNela (see Table 4). Here, 
in the case of NGC 1569, we discuss results from models in which 
both 77SNH and 7/sNia are changed. In two cases, also the prescrip- 
tions about the DM content are modified. This affects the depth of 
the potential well, and thus the time of the onset of the galactic 
outflow. 

In Table 7 we list results from models analogous to Model 5aN 
(Table 7, first row) as far as the prescriptions about the SFH, stel- 
lar nucleosynthesis and GW efficiencies ft are concerned, but in 
which different values for t/snii and 7/sNia are adopted. In the ex- 
treme case where t/snii is set to 0.03 and r/sNia is set to 1, the time 
of occurrence of the GW is delayed (tew — 4 Gyr rather than 
2.5 Gyr; Table 7, last col, first vs. second row). The net outflow rate 
is reduced, since more weight is given to the less numerous SNela, 
according to Equation (7). As a result, we predict a larger present- 
day gaseous mass, a higher ISM metallicity at the present time, and 
a lower mass loss rate, which does not track the SFR any more. This 
is due to the fact that, while the SNII rate closely follows the SFR, 
the SNIa rate displays a very different trend. Therefore, while the 
energy injection from SNell proceeds in lockstep with the SFR, the 
bulk of the explosion energy from SNela is restored to the ISM with 
a delay time which depends on the properties of the SF burst. For 
Model 5aJV, this delay time is ~ 2 Gyr for the early weaker activ- 
ity and ~ 1 Gyr for the second to last burst. SNIa progenitors born 
during the most recent SF activity had no time to restore the bulk 
of their energy into the ISM yet. If one increases the values of the 
ft parameters - which are free parameters of the model - in order 
to reproduce the observational constraints, finds that by assuming 
higher ft values (~5x 10~ 6 for H, D and He and ~ 5 x 10~ 5 
for the elements heavier than helium) can fit all the available ob- 
servational constraints, but the present-day mass loss rate, which 
turns out to be roughly a factor of three lower than the observed 
present-day SFR. This seems to suggest that 7/snii — J/SNia is a 
more reasonable choice. 

In Table 7 we also show the results obtained in a case where 
t/snii = »?SNia = 0.50. In this case, a lower present-day gaseous 
mass survives the SF process (fourth col., third row), unless one 
reduces the strength of the ft coefficients. It is very difficult to put 
constraints on the values of the ft and jjsnii, ??SNia parameters sep- 
arately. In fact, they enter the model in a complex way, i.e. through 
the products ft x t/snii, ft x r/sNia, which are then used to weigh 
the SN rates. The poorly known DM amount and distribution do 
not seem to cause further complications: in Table 7 we display the 



results obtained by a model in which Md/Mi um is set to 50 rather 
than 10 and t/snii = T/SNia = 0.50 (fourth row). It can be seen that, 
although in this case £gw — 2 rather than ~ 1 Gyr, the model pre- 
dictions are almost the same (cfr. third vs. fourth rows in Table 7), 
due to the scarce evolution of the system at those early times. Sim- 
ilarly, setting M d /Mi um to 5 and t/snii = T/sNia = 0.20 does not 
affect the model predictions (Table 7, last row). 

4.2.3 Late accretion of gaseous lumps onto NGC 1569? 

It has been suggested (Melioli & de Gouveia Dal Pino 2004) that in- 
fall of fresh gas from the surroundings could help to keep the value 
of the thermalization efficiencies lower than unity. Therefore, we 
have computed a model (Model 5dN) in which the infall rate in- 
creases with time and the time scale for infall, r, is set to 5 rather 
than 0.5 Gyr, simulating a late accretion of gaseous lumps. With 
this choice, the present-day properties of NGC 1569 are well repro- 
duced (see Fig. 8). The unusually high SFR in NGC 1569 is now 
explained by the replenishment of the HI reservoir in its disc by 
newly accreted material. Indeed, some evidence that the galaxy is 
presently ingesting a small gaseous companion exists (Stil & Is- 
rael 1998; Miihle et al. 2005). In particular, according to Miihle et 
al., the emergence of the latest burst could be explained by the ac- 
cretion of HI gas with a mass of ~ 10 7 M Q . Our Model 5dN is 
characterized by an infall rate Mint ~ 0. 1 M yr _1 during the past 
100 Myr, which implies that during this period nearly 10 7 M© of 
primordial gas are accreted by the galaxy. The predicted outflow 
rate is Mgw ~ 0.6 M Q yr _1 , which leads to a total ejected mass 
of 6.3 x 10 6 Mq over the last 10 Myr, to be compared with an up- 
per limit to the wind mass of 6.2 x 10 6 M obtained from spectral 
fitting (Martin et al. 2002). 

4.3 The problem of nitrogen: NGC 1569 vs NGC 1705 

Understanding the evolution of nitrogen is one of the major chal- 
lenges of modern astrophysics. The situation is far from being 
settled even in our own Galaxy. Recent abundance data for very 
metal-poor Galactic halo stars (Israelian et al. 2004; Spite et al. 
2005) suggest that an important production of primary N took place 
in the massive halo stars, while the delayed N production from 
intermediate-mass stars likely overwhelms any massive star contri- 
bution at later times (e.g. Chiappini, Romano & Matteucci 2003a). 
Yet, current massive star models do not produce enough nitrogen 
to fit the halo data (see Chiappini, Matteucci & Ballero 2005 for an 
extensive discussion). 

Here we find that, while the N/O ratio in the local ISM of 
NGC 1569 is rather well reproduced by models assuming ftj/ft) = 
1, in order to explain the low N/O ratio measured in HII regions 
belonging to NGC 1705, it is necessary to assume ftj/ftj ^ 1. 
This result is independent of the adopted stellar yields and infall 
law. Furthermore, when steepening the IMF slope in NGC 1705 
as suggested by Annibali et al. (2003; x = 1.6 rather than 1.35), 
the discrepancy between model predictions and observations in- 
creases. Before concluding that GWs in NGC 1705 must eject ni- 
trogen more effectively than in NGC 1569, there is one more pos- 
sibility we have to examine. 

4. 3. 1 IMF variations from galaxy to galaxy ? 

One may wonder whether coupling an IMF different from Salpeter 
in NGC 1569 with the same ftj /ft) value which matches the N/O 
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Figure 8. Predictions of Model 5dN for NGC 1569: a) M gas , b) SFR, c) iW s tars> d) gas metallicity, e) Type II (solid line) and Type la (dashed line; multiplied 
by a factor of 10 to make it clearly visible) SN rate, f) Mqw (solid line) and Af; n f (dashed line) as functions of time; g) log(0/H)+12 and h) log(N/0) as 
functions of time. The stars represent observed present-time values (see text for details). 



data for NGC 1705 would allow to reproduce the N/O ratio mea- 
sured for NGC 1569 as well. Possible IMF variations in the range of 
massive and intermediate-mass stars are constrained by HST pho- 
tometry. In particular, for NGC 1569 a power-law exponent flatter 
than Salpeter's cannot be excluded. Indeed, a value of x = 0.35 
over the whole stellar mass range still reproduces the red plume 
of the observed CMD (Angeretti et al. 2005), though it should 
be cautioned that, due to NICMOS resolution limits, small unre- 
solved stellar clusters might appear as single stars populating the 
red plume, thus leading to an artificial flattening of the IMF. Be- 
sides this, one should always bear in mind that the observed CMDs 
cannot constrain the IMF below the mass limit given by the galaxy 
distance and the photometric depth of the images. 

In the following, we show results obtained by running modi- 
fied versions of Model 5aJV, which assume different IMFs, still full- 
filling the required photometric properties of NGC 1569. We adopt 
either x = 0.35 over the whole stellar mass range, or a two-slope 



IMF consistent with Salpeter down to 0.5 Mq, and then flattening 
to x = 0.3 for 0. 1 < m/M < 0.5 (Kroupa & Weidner 2005). 

The input SFRs are modified accordingly (a straightforward 
scaling law does not apply - see Angeretti et al. 2005, their fig- 
ure 6). Extrapolating the SFR below 0.5 M© by adopting x = 0.3 
rather than x = 1.35, results in a SFR about 30 per cent lower than 
assumed up to now (Angeretti et al. 2005). As a consequence, a 
lower final mass in stars of M s tars(tnow) ~ 8 x 10 7 Mq and a 
higher present-day gaseous mass of A/ gas (t now ) ~ 2.5 x 10 s Mq 
are predicted by Model 5aJV. The predicted current metallicity of 
the galaxy is Z = 0.003, while the theoretical N/O ratio is almost 
unchanged, log(N/0) = —1.33 at the present time (cfr. Fig. 5, last 
col, last row to Fig. 9, bottom panel, dashed line). However, since 
the SF preceding the last 1-2 Gyr activity unraveled by HST is a 
poorly-constrained parameter of the model, we can enhance it with 
respect to the prescriptions of Model 5aN. For example, by set- 
ting v = 0.02 Gyr -1 rather than 0.012 Gyr -1 , we find a final stellar 
(gaseous) mass of 1 x 10 8 M (1.8 x 10 s M ). The final metallic- 
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Figure 9. Temporal evolution of oxygen (top panel) and nitrogen to oxygen 
(bottom panel) in the ISM of NGC 1569. Model predictions for different 
IMF slopes [x = 1.35 (0.30) above (below) 0.5 M , dashed lines; x = 0.35 
for the whole stellar mass range, dotted and solid lines) and metal ejection 
efficiencies (/3; = 2 X 10~ 5 , dashed and dotted lines; ft ~3 X 10 -5 , solid 
lines) are compared to the available observations (filled squares; see text). 
The ending points of the evolutionary tracks are marked with arrows in the 
log(N/Q) vs. time plot. 



ity of the ISM is Z = 0.004. Therefore, a two-slope IMF consistent 
with Salpeter's one above 0.5 Mq and flattening below (x = 0.3) is 
permitted by the models. 

If we assume instead x = 0.35 over the whole stellar mass 
range, we find a stellar (gaseous) mass of 6 x 10 6 Mq (2.9 x 10 8 
Mq), whereas a metallicity as high as Z = 0.006 is attained, so 
the model cannot be accepted. Notice that by increasing the metal 
ejection efficiency (i.e., the f3{ parameter value) one can recover the 
present-day metallicity, but the predicted nitrogen to oxygen abun- 
dance ratio is always out of the range allowed by the observations. 
This can be clearly seen from Fig. 9, where the dotted lines stand 
for a model analogous to Model 5aN, in which the IMF slope is 
set to 0.35 rather than 1.35, while the solid lines refer to the same 
model once the metal ejection efficiency is increased in such a way 
that Z = 0.004 at the present time. 

In conclusion: (i) IMF variations compatible with the shape 
of the CMD do not prove useful in solving the puzzle of the dif- 
ferent log(N/0) measured for NGC 1569 and NGC 1705; (ii) flat- 
tening the IMF for 0.1 < m/M© < 0.5 in NGC 1569 produces 
results which still agree with the observations, whereas (iii) adopt- 
ing a flat IMF over the whole stellar mass range produces results at 
variance with the observational data. This may be a hint that small 
unresolved stellar clusters have really been treated as single stars 
located on the red plume in the CMD analysis. We thus conclude 
that different efficiencies of nitrogen ejection through the winds are 
needed in order to explain both the NGC 1569 and NGC 1 705 rele- 
vant data. 



5 DISCUSSION AND FUTURE WORK 

In this section we extend the previous discussion to the interpreta- 
tion of the log(N/0) vs. log(0/H)+12 diagram for a wider sample 
of DIGs and BCDs. In Fig. 10, upper panel, we show galaxies for 
which both oxygen and nitrogen have been measured. The posi- 
tion of NGC 1705 and NGC 1569 in the diagram is given by the 
big filled diamonds^ and the big filled square with 2-cr error bars, 
respectively (see Sect. 2.2 for references). We remind the reader 
that this is not an evolutionary diagram: it refers to present-time 
abundances, derived from H II region spectroscopy. 

Several features emerge from this diagram. It has been pointed 
out (Izotov & Thuan 1999) that the N/O ratio is independent of 
oxygen for BCDs with log(0/H)+12 < 7.6. As already discussed 
by Matteucci & Tosi (1985), this is easily interpreted if the frac- 
tion of primary N is at least 30 per cent. Izotov & Thuan interpret 
this result as an indication for primary nitrogen being produced by 
the same massive stars which synthesize oxygen. In their view, the 
small dispersion in N/O at these low metallicities argues against 
any delayed primary nitrogen production from intermediate-mass 
stars, thus suggesting that these galaxies are now undergoing their 
first burst of SF. The delayed release of primary N would not leave 
its imprints until metallicities 7.6 < log(0/H)+12 < 8.2 are at- 
tained, while secondary N production would become important 
only when galaxies get enriched to log(0/H)+12 > 8.2 (see also 
Henry, Edmunds & Koppen 2000). These mechanisms would be 
responsible for the scatter in N/O observed at log(0/H)+12 > 7.6. 

In previous works (Matteucci & Tosi 1985; Marconi et al. 
1994; Chiappini et al. 2003a) we were able to populate the observed 
log(0/H)+12-log(N/0) diagram with model galaxies, without in- 
voking a primary N production from massive stars (see also Henry 
et al. 2000). We found that dwarf galaxies experiencing an initial 
burst of SF in the past and a more recent one, or, alternatively, just 
one SF episode lasting long enough to let intermediate-mass stars 
have time restoring their nucleosynthesis products, tend to cluster 
around the N/O 'plateau'. Models in which more than one burst 
occur populate the high metallicity, high N/O part of the plot (e.g., 
Chiappini et al. 2003a, their figure 13 and table 5). However, in 
those models the SFH and IMF were treated as free parameters, a 
simpler dependence of the outflow rate on the SFR was assumed, 
and no attempts were made at modelling specific objects. 

As far as an important production of primary N from very 
metal-poor massive stars is concerned, we recall that current mas- 
sive star models do not predict an amount of primary N sufficient 
to explain the Galactic halo star data (Chiappini et al. 2005, and 
references therein). On the other hand, the amount of primary N 
synthesized in low-metallicity intermediate-mass stars due to the 
third dredge-up and HBB processes has been probably overesti- 
mated (e.g. Chiappini et al. 2003a). Besides the observational data, 
in Fig. 10 we plot also the predictions from Models JaM and lfH 
for NGC 1705, and 5aM and 5dH for NGC 1569. In the upper 
panel, only the ending points of the theoretical tracks are shown 
(big crosses), whereas in the lower panel the complete evolutive 
tracks are displayed. NGC 1705 and NGC 1569 show very high 
SFRs - among the highest ever observed for DIGs and BCDs - 
and evidence for protracted SFHs in the past. Therefore, it is not 
surprising that they reach a relatively large metallicity. However, 



^ For NGC 1705, both the mean N/O from measurements in all H n regions 
and that without the extremely low nitrogen abundance of region B4 are 
considered (see Lee & Skillman 2004). 



© 2005 RAS, MNRAS 000, 1-20 



18 D. Romano et al. 



-O.E 



-1.2 

O 

^ 1 A 

O —1.4 



-1.6 

-1.8 
-1.2 




NGC 1569 



5dH 



Jffl 



NGC 1 705 



o 



-1.4 
-1.6 
-1.8 




Figure 10. Upper panel: log(N/0) vs. log(0/H)+12 diagram for several 
samples of DIGs, BCDs, and low surface brightness dwarf galaxies. Filled 
diamonds: metal-poor Hii galaxies from Kobulnicky & Skillman (1996); 
filled circles: low surface brightness dwarfs from van Zee et al. (1997): big- 
ger circles are objects with higher average SFR in the past; the big empty 
triangle stands for the sample object with the highest SFR at present — 0.35 
M0 yr _ 1 , while the small empty one stands for the object with the low- 
est SFR at present - 0.0015 Mq yr _1 ; stars: ground-based spectroscopic 
observations of 54 supergiant H n regions in 50 low-metallicity BCDs from 
Izotov & Thuan (1999); squares: Virgo DIGs and BCDs from Lee et al. 
(2003b; small filled squares) and Vflchez & Iglesias-Paramo (2003; small 
empty squares). The big diamonds and the big square on the right represent 
measurements for NGC 1705 and NGC 1569, respectively, along with their 
2-<j errors. Also shown are predictions from Models iaM, 5aM, lfH and 
5dH (final points of the tracks; big crosses). A higher efficiency of nitrogen 
ejection through the wind is assumed for Model lfH with respect to Mod- 
els JaM, 5aM and 5dH (see text for details). Lower panel: full evolutionary 
tracks for Models iaM, 5aM, lfH and 5dH in the log(N/0)-log(0/H)+12 
space. 



ses (see previous sections). The case with minimum HBB extent 
well matches the log(N/0) value actually observed for NGC 1569. 
However, in order to make the model predictions agree with the 
observations for NGC 1705 too, one must require that some other 
mechanism(s) is (are) at work. We argue (Sect. 4.1) that differen- 
tial galactic outflows in which different fractions of nitrogen are 
lost with respect to oxygen offer a viable solution. For this rea- 
son, Model lfH has been computed by assuming /3n//9o — 2 
rather than 1, namely, postulating a higher nitrogen ejection effi- 
ciency in the outflow of NGC 1705. Alternatively, in NGC 1705 we 
might be seeing localized self-pollution from newly-born massive 
stars. However, the constancy of the oxygen abundance measured 
in several HII regions of NGC 1705 poses serious challenges to 
this interpretation (see Sect. 4.1.1.2). Using the yields from rotat- 
ing star models by Meynet & Maeder (2002) also allows us to fit 
the data (Fig. 10; Models iaM and 5aM). However, in this case 
only the highest N/O value for NGC 1705 and the lowest one for 
NGC 1569 can be explained, rather than the difference between 
the two. When faced with the whole data set, we think that the 
observed scatter of log(N/0) values at relatively high metallicities 
needs the hypothesis of differentially enriched GWs in order to be 
explained. Clearly, only detailed study of a larger number of in- 
dividual dwarfs, for which the SFH and the IMF had been previ- 
ously inferred by means of CMD analyses, would allow us to un- 
derstand the driving processes which lead to the scatter. Recently, 
our group (PI: A. Aloisi) has been granted time to perform HST ob- 
servations of NGC 4449, a starburst DIG located at a distance D ~ 
4.2 Mpc with log(0/H)+12 = 8.3 and log(N/0) = -1.2 (Sabbadin, 
Ortolani & Bianchini 1984), i.e., with a N/O ratio higher than both 
NGC 1569 and NGC 1705. The use of the Advanced Camera for 
Surveys (ACS) on board HST will allow us to reach stars quite 
fainter than the red giant branch tip (if they exist), and to infer the 
SFH of NGC 4449 up to several Gyrs ago. We want also to call at- 
tention to IIZw40, a galaxy located at a distance D ~ 7.5 Mpc - 
i.e., within reach of HST instrumentation - with log(0/H)+12 = 



8.091] 



log(N/0) = 



-1 052+ 059 

'• UJi -0 .077' 



He/H = 0.093 ±0.013 



(Kobulnicky & Skillman 1996), and an old underlying stellar popu- 
lation (Vanzi et al. 1996). Given its very high N/O abundance ratio, 
it would be an interesting object to study in the context of our re- 
search. 



6 CONCLUSIONS 

In this paper we have presented one-zone chemical evolution mod- 
els for two dwarf starburst galaxies, NGC 1705 and NGC 1569. For 
both objects, we have adopted the SFH and IMF previously inferred 
from HST optical and NIR CMDs and found that strong galactic 
winds develop owing to their violent SF activity. Our main conclu- 
sions can be summarized as follows: 



NGC 1705 has the lowest N/O ratio observed at those metallicities, 
and the N/O ratio of NGC 1569 is quite low too, in contrast with 
previous ideas that DIGs and BCDs evolved through several bursts 
of SF or experiencing a long-lasting SF regime with a high SF ef- 
ficiency should populate the high N/O region of the log(0/H)+12- 
log(N/0) plot. Models lfH and 5dH have been computed with the 
van den Hoek & Groenewegen (1997) plus Nomoto et al. (1997) 
yields, for the case of minimum HBB extent. Both assume an 
infall law exponentially increasing in time. The adopted stellar 
IMF (Salpeter 1955) and SFHs are consistent with CMD analy- 



(i) The current metallicity of HII regions of both NGC 1705 
and NGC 1569 is well reproduced only if strong GWs efficiently 
remove the metals from the galaxies, in agreement with previous 
works (e.g., Matteucci & Tosi 1985; Pilyugin 1993; Carigi et al. 
1995; Recchi et al. 2001, 2004, 2005). 

(ii) The very low N/O ratio measured in NGC 1705 (Lee & 
Skillman 2004) requires a higher efficiency of nitrogen ejection in 
the wind of NGC 1705 relative to the wind of NGC 1569 in or- 
der to be explained, a result independent of the other parameters 
of the model. Alternative solutions, such as complete suppression 
of HBB in intermediate-mass stars or localized chemical pollution 
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from very massive stars in NGC 1705, are discussed and seem un- 
likely. 

(iii) An offset exists between the oxygen abundance of the neu- 
tral phase and that of the ionized gas in NGC 1705, which seems 
absent in the N/O data (Aloisi et al. 2005). Once some primary 
nitrogen production through HBB in intermediate-mass stars is al- 
lowed for in NGC 1705, an almost constant N/O abundance ratio 
is predicted in its ISM during almost all the evolution, whereas the 
oxygen abundance always increases with time (see, e.g., Fig. 3, first 
to third cols.). This qualitatively explains why an observational off- 
set should be present in O/H, and absent in N/O, though uncertain- 
ties in both the physics of the models and the treatment of the data 
prevent us from attempting a more quantitative analysis. Notice that 
by using the yields of Meynet & Maeder (2002), computed by tak- 
ing into account mass loss and stellar rotation, but neglecting HBB, 
the theoretical N/O ratio in NGC 1705 steeply increases during its 
whole evolution (Fig. 3, last col., last row). Therefore, in this case 
the constancy of the N/O ratio in the neutral and ionized gas phases 
cannot be explained. 

(iv) We find hints that the third dredge-up and HBB must be less 
effective than assumed up to now in low-metallicity intermediate- 
mass stars. We do not find any need for a significant primary nitro- 
gen production from massive stars to explain the present-day abun- 
dances in NGC 1705 and NGC 1569. The large primary nitrogen 
production which characterizes massive stars at Zi n i < 10~ 5 (see 
Chiappini et al. 2005) does not leave any imprint on the present-day 
properties of NGC 1705 and NGC 1569, at Z ~ 0.004. 

(v) We demonstrate that the efficiency of hydrogen entrainment 
in the outflow crucially depends on the adopted infall law. If we al- 
low for important mergers with gaseous lumps at late times, we find 
that the efficiency of hydrogen ejection through the outflow must 
be of the same order of magnitude as that for the heavy elements 
in order to fit the observations. On the contrary, with an infall rate 
exponentially decreasing with time (as usually assumed by GCE 
models), the required efficiency of hydrogen ejection is up to one 
order of magnitude lower than that for the heavy elements. While 
some evidence exists that NGC 1569 could be currently accreting a 
small companion cloud (Stil & Israel 1998; Miihle et al. 2005), no 
clear signs of external triggers to the sudden, intense ongoing burst 
observed in NGC 1705 have been found up to now (Meurer et al. 
1998). 

(vi) By assuming that the GW efficiency varies in lockstep with 
the predicted Type II plus Type la SN rate, we find Mgw (t now ) ~ 
ip(t now ) in NGC 1569 [where MGw(inow) and ip(t now ) are the 
current mass loss and SF rates], independently of the assumed in- 
fall law and in agreement with observations. The total mass of 
gas ejected at the average rate predited during the last 10 Myr is 
6.3 x 10 6 M in the case of an infall rate exponentially increas- 
ing with time, while it is 3.5 x 10 6 Mq in the case of an expo- 
nentially decreasing infall. These values agree with the best-fitting 
value (3.5 x 10 6 M Q ) and upper limit (6.2 x 10 6 M ) for the 
mass ejected in the outflow given by Martin et al. (2002). Cur- 
rent estimates of the gas mass in the outflow of NGC 1705 vary 
by almost three orders of magnitude, going from Mgw(£iiow) — 
0.0026 Mq yr _1 up to ~ 5 Mq yr _1 (Meurer et al. 1992). These 
values do not usefully constrain our models, which predict mass 
loss rates ranging from MGw(inow) — 0.04 to 0.2 M yr _1 , de- 
pending on the adopted infall law prescriptions. 

Clearly, we need to apply our models to a larger galaxy sample 
before drawing firmer conclusions about late-type dwarf galaxy 
formation and evolution. We identify NGC 4449 and IIZw40 as 



promising candidates: owing to their metallicities, similar to those 
of NGC 1705 and NGC 1569, but higher N/O ratios, they probe a 
different locus in the log(N/0) vs. log(0/H)+12 plane. Thus, they 
could help us to unravel the origin of the observed spread among 
N/O values at a given metallicity. In particular, our group has al- 
ready been granted HST time to perform ACS observations of 
NGC 4449. These will enable us to determine its SFH and IMF, 
a necessary step to work out the appropriate chemical evolution 
models. 
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